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Abstract 

A novel continuum theory of incoherent interfaces with triple junctions is applied to study 
three-dimensional coupled grain boundary (GB) motion in polycrystalline materials. The kinetic 
relations for grain dynamics, relative sliding and migration of the boundary, and junction evolution 
are developed. In doing so a vectorial form of the geometric coupling factor, which relates the 
tangential motion at the GB to the migration, is also obtained. Diffusion along the GBs and 
the junctions is allowed so as to prevent nucleation of voids and overlapping of material near the 
GBs. The coupled dynamics has been studied in detail for two bicrystalline and one tricrystalline 
arrangements. The first bicrystal consists of two rectangular grains separated by a GB, while the 
second is composed of a spherical grain embedded inside a larger grain. The tricrystal has an 
arbitrary shaped grain embedded inside a much larger bicrystal made of two rectangular grains. 
In all these cases, analytical solutions are obtained wherever possible while emphasizing the role of 
various kinetic coefficients during the coupled motion. 

Keywords: Incoherent interfaces; Triple junction; Coupled grain boundary motion; Geometric 
coupling factor; Nanocrystalline materials 


1 Introduction 

We develop a thermodynamically consistent continuum framework to study three-dimensional (3D) 
coupled grain boundary (GB) motion in the presence of triple junctions. A GB is modelled as a sharp 
incoherent interface connected to other GBs at junction curves. The irreversible dynamics at a GB is 
governed by its normal motion (GB migration) and a relative tangential sliding of the adjacent grains. 
The latter can arise due to the inter-granular viscous sliding, possibly as a result of the twist component 
of the GB, or/and as a result of coupling with GB migration [S]. In polycrystalline materials with rigidly 
deforming grains, as will be assumed presently, the sliding can be decomposed into a relative translation 
and a relative rotation between the adjacent grains. On the other hand, the irreversible dynamics at 
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a junction is governed by the motion of the non-splitting junction curve. The presence of junctions 
can significantly influence the overall dynamics of all the GBs and the grains in their neighborhood, 
for instance by inducing drag or altering diffusive flux [T]. The coupled motion, which requires sliding 
to be necessarily coupled with GB migration, is the dominant mechanism for both grain coarsening 
and plastic deformation in nanocrystalline (NC) materials with average grain size of the order of few 
tens of nanometers (hence a large volume fraction of GBs and triple junctions) [16[l23j . This is unlike 
coarse-grained materials where GB migration and dislocation dynamics dominate grain coarsening and 
plastic deformation, respectively. The coupled motion has recently been studied theoretically [Tll2l[5lin] , 
experimentally m, and with molecular simulations [22] . Although some of these studies have included 
the effect of junction dynamics [I1E2], all of them are restricted to two-dimensional grains and hence 
applicable only to polycrystals where each grain is columnar and identical in cross-section along the 
length direction; such a restriction requires the GB to have only tilt, and no twist, character. 

The main contributions of this paper include: 

(i) A 3D thermodynamic formalism including diffusion to deal with incoherent interfaces with junctions 
(Section [3|). Junctions have been previously studied in the context of continuum thermodynamics but 
only with coherent interfaces and without diffusion |6l[20]. On the other hand, thermodynamics of 
incoherent interfaces has been explored earlier without considering junctions |7|- All of these works 
were based on the framework of configurational mechanics. Our treatment, while extending to junctions 
with incoherent interfaces, takes an alternate viewpoint where we do not regard the configurational 
forces to be fundamentally on the same footing as standard forces (with their own balance laws etc.). 
We introduce configurational forces in our formalism as mechanisms of internal power generation so as 
to ensure that the excess entropy production is restricted to interfaces and junctions. A 2D version of 
this formalism was recently presented by the authors [T]. 

(ii) Deriving kinetic relations for coupled GB motion in three dimensions (Section [5|). We extend earlier 
models of coupled GB motion to a 3D setting. The first kinetic relations for the coupled motion were 
proposed by Gahn and Taylor nm which were restricted to two-dimensions and only bicrystalline 
arrangements (hence no junctions). They also ignored the possibility of relative translation of grains 
while considering sliding at the GB only due to the relative rotation. More recently, the present authors 
have extended the model to include junctions and relative translation but still restricting themselves 
to two dimensions pQ. 

(iii) Formulation of a vectorial geometric coupling factor (Section Uj). The coupling between the tangen¬ 
tial and the normal motion of the GB is purely geometric and depends on the measure of incoherency 
at the boundary |3H5]. The incoherency is quantified by the net Burgers vector (given by Frank-Bilby 
relation) or equivalently by the interfacial dislocation density. The GBs in the present 3D framework 
generally have a mixed character with both tilt and twist components. The coupling factor for a high 
angle planar symmetric tilt boundary, derived previously by Gahn et al. jam, therefore needs to be 
extended to include multiple sets of edge and screw dislocation arrays. The coupling factor now derived 
is a vectorial quantity rather than a scalar as has been the case in the earlier studies. 

Our derivation for kinetic relations is based on the following assumptions: (a) the individual grains 
experience only rigid deformations (i.e. translations and rotations), (b) the shape accommodation 
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required for preventing void-formation and interpenetration of the material in proximity of the GBs, 
during relative tangential motion between the grains, is accomplished by diffusion across as well as 
along the GBs and also along the junction curves; (c) the velocities associated with various GBs, 
grains, and junctions remain much smaller than the speed of sound in the material. The inertial 
effects are therefore ignored; (d) the grains are free of defects and all the lattice imperfections are 
concentrated at the GBs and junctions (this is reasonable for NC materials with their small grain size); 
and (e) no additional stress fields are present at the interface and the junction. The GBs are considered 
to be orientable surfaces (of arbitrary shapes) with five macroscopic degrees of freedom which include 
three misorientation angles and two independent variables describing the orientation of the GB. The 
junctions are arbitrary 3D space curves with varying curvature, normal, binormal, torsion etc. The 
excess energy density of a GB is assumed to depend on the five parameters mentioned above, while the 
excess energy density of a junction is assumed to depend only on the unit tangent associated with the 
junction curve. 

The paper has been organized in the following manner. In Section [2] we briefly introduce various 
kinematic and integral relations required for our study. We derive the essential balance laws and 
dissipation inequalities in Section [3l A generalized derivation of the vectorial geometric coupling factor 
has been presented in Section 01 In Section [5] we apply our theory to derive the kinetic relations for GB 
motion, grain dynamics, and junction motion for two bicrystalline and one tricrystalline arrangements. 
The phenomenological kinetic equations are motivated from the dissipation inequalities derived from 
the second law of thermodynamics in confirmation with other standard balance laws of continuum 
physics. One bicrystal has two rectangular grains separated by a low angle planar mixed GB, while the 
other has a spherical grain embedded inside a larger grain. The tricrystal constitutes of an arbitrary 
3D grain which is completely embedded inside a bicrystal consisting of two large rectangular grains. 
We conclude our study with a discussion on some open directions in Section [ 6 l 


2 Kinematics 


We consider a region P, as shown in Figure [l(b)[ taken out from a polycrystalline arrangement depicted 
in Figure 1(a) It contains three subgrains Pi, P 2 , and P 3 , three smooth GBs Fi, r 2 , and Fa, and a 
smooth junction curve J. The normal n* to Fj is chosen such that it points inside Pi, where i = 1, 2, 3. 
We denote the position vector of a point by x and the time by t. The grains are oriented differently 
with respect to a fixed coordinate. 

Let A and B be two second order tensors. The derivative of a scalar valued differentiable function 
of tensors, say G{A), is a tensor BaG defined by 


G{A + B) = G{A) + dAG-B + o{\B\), (1) 

where o{\B\)/\B\ —>■ 0 as \B\ —>■ 0; the norm of a tensor is defined as |Sp = B ■ B. Similar definitions 
can be made for vector and tensor valued differentiable functions (of scalars, vectors, and tensors). 
The derivative of a field defined over P with respect to the position vector is denoted by the gradient 
operator V. 
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(a) 


Figure 1: (a) Schematic of a polycrystal in 3D. (b) The region P containing three subregions Pi, P 2 , 
and P 3 . The GBs OO'A!OO'B'B, and OO'C'C (with normals rii, n 2 , and as shown), are 
denoted by Ti, r 2 , and Ts, respectively. The curve 00' is the triple junction J. 

2.1 Bulk fields 


Let / be a piecewise smooth bulk field which is discontinuous across Tj and singular at J. We denote 


the jump of / across Tj as |/] = / , where /+ is the limiting value of / at a; £ Tj from the 

grain into which rij points and f~ is the limiting value from the other grain. If /i and /2 are two 
piecewise continuous functions across Tj, then I/1/2I = [/il(/2) + (/i)[/2l, where (/) = (/+ + f~)/2 
is the average of /'*' and f~. To deal with the singularity of the field at the junction we carry out 
our analysis in a punctured region P^ obtained by excluding a small tube Tg of radius e from P in the 
neighborhood of the junction, cf. [5U]. The outward normal to the boundary of the tube dT^ is denoted 
by m. The boundary dT^ moves with a velocity u. 

We assume that / satisfies the limit 



( 2 ) 


where dv is an infinitesimal volume element of P. Using the standard transport relations for the bulk 
quantities it can be shown that [ 20 ] 



where the superposed dot denotes the material time derivative, div is the divergence operator, v is the 
particle velocity, Vi is the interfacial normal velocity, Ui = Vi-v-rii is the relative normal velocity of the 
interface, and da is an infinitesimal area element of a surface. Let a and A denote piecewise-smooth 
vector and tensor fields, defined in P, which are singular at the junction. The divergence theorem 
requires m 



( 4 ) 
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diY Adv = / Amda — ^y / \A\nida — \im / Amda. 

' JdP Jvi <^^0 JdT^ 


(5) 


2.2 Interfacial fields 


Consider an orientable surface T (subscript i is presently dropped), with boundary dF, and let n and 
V be the associated unit normal field and normal velocity field, respectively. The surface gradient of a 
scalar field g, vector field g, and tensor field G, all smoothly defined over F, are defined as 

= P{Vg), = {Vg)P, and = (VG)P, ( 6 ) 

respectively, where P = I — ni^n is the projection tensor (/ is the 3D identity tensor and 0 denotes 
the dyadic product); while calculating Vg (etc.) one has to use a smooth extension of 5 in a small 
neighborhood of F. The surface divergence of these helds are defined by 

diV^^ = tr(V^g) and k ■ div^ G = div^(G'^fc), (7) 

for all constant vectors k, where tr represents the trace operator and the superscript T stands for the 
transpose. The surface Laplacian of g is given by 

A^g = div^iV^g). ( 8 ) 

The curvature tensor field L and the total curvature k associated with F are defined as 

L = — V^n and k = tr L, (9) 


respectively. 

Let t be the outward unit normal to the closed curve clF such that n-t = 0. When G and g satisfy 
Gn = 0 and g ■ n = 0, respectively, the surface divergence theorem yields m 


/ g ■ tdl = di 

Jdr Jr 


/ Gtdl = / div'^ Gda and 

Jdr Jr 

where dl is an infinitesimal line element along a curve. 

The normal time derivative of g following F is given by m 


div* gda, 


( 10 ) 


g = g + VVg ■ n, 


( 11 ) 


which is the rate of change of g as experienced by an observer sitting on the moving surface F. The 
first term indicates the local rate of change of g at a fixed material position, while the second term 
represents the rate of change of g due to influx of particles along n as the interface moves with velocity 
V. The following identities can be readily verified m- 

n = -V^V andL = -V^n-Ln0n + VL^. (12) 

On the other hand, the intrinsic time derivative of g following (9F is given by [U] 

g = g + Vg -w, (13) 
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where w is the intrinsic velocity of 9r, such that 


w = Vn + Wt, (14) 

and W is the velocity of clF along t. According to (I13p the rate of change of g following OF is equal 
to the sum of the rate of change of g following F and a term representing the change in g due to the 
incoming particles from the neighborhood along the tangential direction t. 

We will need the following transport theorem for F such that a part of 9F intersects with dP and 
the rest with J m- 

^ f gda= [ (g - gKV)da+ [ gWdl+ [ gq -tdl, (15) 

dt Jr Jr JrndP Jj 

where is the intrinsic (independent of the parametrization) velocity of the junction. 


2.3 Junction fields 


Let 6 and b denote the the terminal points of the junction curve J, and let I be the unit tangent to 
the curve such that it is directed towards b. The normal and binormal vectors associated with J are 
denoted by u and b, respectively. The projection tensor Q = I — — b®b = l®l maps any 

vector on to the tangential direction of J. The intrinsic velocity field of the junction, which can be 
decomposed as 

Qp = Quf + Qbb, (16) 

is such that {I — Q)u —as e —)■ 0 [2U]. The velocity of terminal points is denoted as q such that 

Q = qp + ¥ (17) 


at the respective end points. The intrinsic time derivative of a scalar field defined on the junction 
curve, say x is given by (compare with (fT3]) i 


X = X +Vx • qp. 


The transport theorem associated with x is given by 

d 


dt Jj ^ ~ ^ (X®)5, 


(18) 


(19) 


where kj is the curvature of the junction curve. The gradient of x along the junction curve is defined 
as V'^x = Q^X- Similarly for a vector field defined on J we introduce V'^qp = {'Vqp)Q. We note the 

identity m 

I = {V^qp)l + q^Kjl. ( 20 ) 

It is useful to decompose an integral over the tube surface around the junction as 


lim 

e-i-O 


ada = 


lim 


/9Te JJ 

where dT^, is the envelope of the circles Cg of radius e. 


e-i-O 


adl 


a 


dl, 


( 21 ) 
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3 Balance laws and dissipation 


We now obtain the consequences of balance of mass and momentum, as well as obtain local dissipation 
inequalities in the bulk, at the interface, and at the junction, while restricting to the assumptions 
enlisted in Section [TJ 


3.1 Balance of mass 

The rate of change of total mass in P is balanced by the mass transport into the region via bulk 
diffusion across dP, GB diffusion at the edge Tj n dP, and diffusion at i and 5. Neglecting excess mass 
densities of the GBs and the junction, the mass balance can be written as 

^ [ pdv = - [ j ■ mda - V' / hi- Udl - {hj)s, (22) 

Jp JdP ^ JVindP 

where p is the mass density of the bulk grain, j is the bulk diffusional flux, hi is the tangential 
diffusional flux on Tj, and hj is the diffusional flux along the junction. Using transport theorem ([3]) 
and divergence theorems (j3|) and m2, and localizing the result owing to the arbitrariness of P, we 
can obtain the following local equations: 

p + pdivu + divji = 0 V® G Pj, (23) 

IpUij = IjI • rii + div*^ hi\fx eVi, and (24) 

lim / p{u — v) ■ mdl = V'^hj • I + lim / j ■ mdl — (hi ■ U) r Vx G J, (25) 

where we have used the limit lim^^o far nr ~ Ij j dl- 


3.2 Balance of linear momentnm 


Neglecting inertia and body forces, and assuming absence of interfacial and junction stress fields, the 
balance of linear momentum is given by 


amda = 0, 


(26) 


JdP 

where a is the symmetric Gauchy stress tensor. Using ([5]) the following local equations are readily 
obtained pO] : 


diver = 0 V® G Pi, (27) 

|cr]ni = 0 Va; G Tj, and (28) 

lim / crmdl = 0 V® G J. (29) 


According to (|28l) the traction field is continuous across the GBs, whereas (1291) requires that the net 
force acting at each circular region is zero in the limit e —>■ 0 , although the stress field can still be 
singular at the junction (weak singularity). 
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3.3 Dissipation inequality 


Let \I' be the bulk free energy density, 7 * the interfacial free energy per unit area of Fj, and rj the free 
energy per unit length of J. For an isothermal environment, the mechanical version of the second law 
requires the rate of change of the total free energy P to be less than or equal to the total power input 
into P PEI], i.e. 


A 

dt 


3 » 

^dv + / 7j(ia + 



< / am ■ vda — fij ■ mda — 

Jap JdP Ja 


apnTi 


fihi ■ tidl 


n 


-{nhj)s + y] / (ci-Wi + n- ni)dl + {u: ■ q)g 
JaPHTi 


(30) 


where /r is the chemical potential. We assume that the chemical potential is continuous across the 
interface and at the junction (i.e. local chemical equilibrium P). The first integral on the right hand 
side of (|30l) is the power input through the tractions acting on dP] the next three terms are contribution 
to power input due to mass flux at dP, dP H Fj, and the end points of J, respectively. The last two 
terms are non-standard; we will discuss their significance before deriving the consequences of (I30p . 
These terms are required to ensure that there is no excess entropy production at the edges dPDTi and 
at the terminal points of J. The excess entropy generation is necessarily restricted to the GBs and the 
junction. These additional power input terms are to be considered in (I30p only when the edges and 
the terminal points lie on the surface of an interior part of a body. The precise form of Cj, Tj, and 
u will depend on the constitutive prescriptions for free energies and stress. At this point these are to 
be understood as agents of power input, in conjugation with the respective intrinsic velocities, so as 
to ensure that the net entropy generation meets the above mentioned requirement. Such terms also 
appear in the framework of configurational mechanics [7l ll4l[20] . where the existence of Cj, r*, and u: 
is assumed a priori as fundamental forces which satisfy certain balance relations. Our treatment (see 
also m^m) is motivated purely from the viewpoint of quantifying excess entropy generation. For 
this we do not have to consider any additional balance laws other than those which are standard in 
continuum physics. 

Using transport theorems ([3|), (E^ . and (fTop . divergence theorems ([1]) and (fTOl) . and the decompo¬ 
sition dm), we can rewrite the inequality pOP as 

6 

^/a<0, (31) 

a=l 


where 


Ii = y div V + pfi — a ■ Vv + j ■ Vpj dv, 

3 . 

h = (li- liKiVi - {UiEjui ■ Ui - {arii) ■ Pi|u] -F hi ■ da, 

i=i 

3 . 

L 

^ JaPHTi 


(32) 

(33) 


— Ci ■ Wi — Ti ■ rii 


dl 


(34) 


h = - 


h — 


/ ( lim / (av + ^{u — v) — p,j) ■ mdl ) dl, 

JjV^^Jc, ) 

J (^- di^JQu + V-\^ihj) ■ I + + 7i9 • dl, and 


h = {mi - ^ 


(35) 

(36) 

(37) 


In obtaining (f32]l we have used (f23]i and (f?7ll . and introduced Vl = ^ — (the grand canonical 
potential). To derive (1331) . on the other hand, we have used (IMl) and (l28l) : here E = QI — <t is the 
bulk Eshelby tensor. 

To determine the precise form of Cj, Tj, and uj, and also to obtain the local dissipation inequalities 
associated with the grains, GBs, and junction, we will now prescribe the constitutive nature of the GB 
energy and the junction energy. Towards this end, we assume the GB energy to depend on the misori- 
entation between the grains, the normal to the GB, and curvature. The former two dependencies are 
standard in material science literature (cf. Chapter 12 in [H]). The curvature dependence is primarily 
introduced to regularize the governing partial differential equations for capillary driven GB motion, 
which otherwise become backward parabolic and hence unstable in certain ranges (GB spinodals) of 
the orientations. We follow Gurtin and Jabbour [T3] in assuming the following quadratic dependence 
of GB energy on curvature: 

7 = 7(0, n) + (38) 

where © is the misorientation tensor given by (i?'*')^i?~; ei and 62 are scalar constants such that 
ei > 0 and €2 + ei/2 > 0. The rotation is the orientation tensor of the grain into which n points, 
and R~ is the orientation tensor of the other grain. We introduce 


M = (9^7 = eiL + €2kP. 


(39) 


It is easy to see that M is symmetric and satisfies MP = M. The junction energy density, on the 
other hand, is assumed to be a function of the unit tangent along the junction curve [20j : 

V = Vil)- (40) 

Substituting (l38]) and (HO]) into (ISTI) . and performing a cumbersome but straightforward calculation, 
yields 

[ Dbdv + / T)rida + [ Djdl - V' / | Wj (7* - MiLtU ■ U - Ci ■ U) 

Jp ^ JTi JJ ^ JdPnVi 

-Vi ((div"^ Mi + dm^i) -ti + Ci- Tii) - Ui ■ {n + Miti)'^dl 

- {mid - m) + {{I - Q)difi - Up) • 9^)5 > 0 , (41) 

where Spi) and T>j are the the entropy generation rates per unit volume of the bulk, per unit area 

of Tj, and per unit length of J, respectively. The expressions for these rates are given in Equations (|46p - 

(1481) below. In deriving the above inequality we have also assumed the junctions to be non-splitting, 
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i.e. Vi = Qp - rii. The first term in the above inequality is the net entropy generation within the grains. 
The next two terms are excess entropy generation at the GBs and at the junction, respectively. The 
rest of terms in (j41j) are the entropy production rate at the edges T* n dP and the terminal points 5 
and i. We however require that the excess entropy production must not have any contribution from the 
edges of the GBs and the terminal points of J, all of which are a part of dP. This is reasonable since 
the entropy generation in P should only be within the grains, at the interfaces, and at the junction. 
Any additional source should vanish. Gonsequently 



Ci ■ ti — 'ji MiLiti ■ ti, 

(42) 


Ci ni = - (diV^ Mi + ■ U, 

(43) 


Ti = —MiU, and 

(44) 


u: = r]l + {I - Q)dirj. 

(45) 

Substituting (I42])-(I45]) back into (jlTTl. and localizing the result, we obtain the following local dissi¬ 
pation inequalities: 


T>b = cr ■ Vu — p + Ddivu) — pfi — j ■ V// > 0 V® G Pi, 

(46) 

Dr, = piEjn 

i-ni + {arii) ■ + fi^i “ (d&P) • ©* > 0 V® € T^, and 

(47) 

Sj = 

Ej-qp — lim / Em ■ vdl — hjiV^p) ■ 1 — Tj • n* > 0 V® € J, 

(48) 

where 

fi = liK-i - div‘^((9ji.7i) - Mi ■ Lf - div‘^(div'^ Mi), 

(49) 


E ] = {I — Q ){lim [ Emdl — c* — / , and 

tt J 

(50) 


f J = -r]Kju + ((J - Q)difi) 1. 

(51) 


Equation (|46l) gives the entropy production rate per unit volume within the grain. For rigidly 
deforming grains cr ■ Vu = 0 and divu = 0. Additionally, if we assume that T = T(/ 3 ) then (I46p yields 
fi = and j ■ V< 0. Equation (1471) contains the dissipation rate per unit area of the GB, with 
contribution from boundary migration, relative translation of grains at the boundary, GB diffusion, 
and misorientation change. The inequality therein forms a basis for postulating kinetic relations for 
coupled GB motion, as is done in Section [5j It can also be a starting point for motivating kinetic 
relations for motion of incoherent phase boundaries with diffusion and curvature dependent boundary 
energy mm, as well as for a variety of physical phenomena involving coherent interfaces [9]. An 
analogous inequality, valid for a one dimensional interface in a 2D grain, was derived recently by the 
authors [ 2 ]. 

Equation (I48p gives the net dissipation rate per unit length of the junction curve, with contribution 
due to motion of the curve, diffusion along it, and evolution of orientation of the intersecting boundaries. 
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A comment is in order regarding the contribution due to the latter, represented by the last term on the 
L.H.S. of the inequality. It is evident from (j39|) and (j44p that this term is linear in scalar parameters 
(ei and 62 ) which appear in the curvature dependent part of the GB energy. The scalar parameters are 
usually infinitesimally small, ensuring that curvature dependent part of the energy is significant only 
at the corners M- The curvatures of intersecting boundaries, as they approach the junction, are finite 
and hence |rj| are small. We can therefore ignore the last term on the L.H.S. of the inequality (|48p 
within the present analysis. Secondly, in the context of GB dynamics, we assume the density field to 
remain bounded at the junction curve and the velocity in each grain to be resulting only a simple rigid 
body motion (hence no strains in the grain). With these assumptions, and keeping in mind the weak 
singularity condition (I29p . we can show that the closed integral terms in (I48|) vanish in the limiting 
sense. We can rewrite (|48l) and (1501) under all these considerations as 

Dj = Tj ■ Qp — hj{V'^fi) • Z > 0 V® G J, where (52) 

:Fj = -{I-Q)(j2c, + fjy (53) 

During thermodynamic equilibrium, with junction curve remaining stationary and diffusion absent, 
= 0 which in the absence of junction energy yields the well known Herring’s relation |15] . i.e. 
iliU ~ = 0. The above framework can be used to obtain the extension of Herring’s 

relation in the presence of junction energy and various singular fields (see also [2D]). Finally, we note 
that junctions have been previously treated in the framework of continuum thermodynamics but only 
for intersecting boundaries which are coherent | 6 l[ 20 ]. The grain boundaries however are in general 
incoherent and a treatment of coupled GB motion necessarily requires allowance for relative slip at the 
boundary. The present framework allows for such incoherency and for junctions which are formed at 
the intersection of such boundaries. 

4 Geometric coupling factor 

During the migration of a tilt or a mixed GB, the adjacent grains undergo a tangential motion giving 
rise to a coupled dynamics (Ghapter 14 in |19j . and [3l ll0l[T7] i. The deformation of the grain in the 
wake of a moving GB, during coupled motion, is essentially controlled by the intrinsic edge dislocation 
content at the GB. Screw dislocations, if present, just glide along the GB plane and contribute only to 
grain sliding without affecting the coupling process [TO]. For a moving planar symmetric tilt GB, whose 
wake experiences a simple shear deformation, Gahn and coworkers [3] introduced geometric coupling 
factor as the ratio of the relative tangential velocity (in the absence of viscous sliding) to the GB 
velocity. For a large misorientation range of a symmetric tilt boundary, containing single array of edge 
dislocations, the coupling factor (denoted by /3) was calculated to be /3 = 2tan(0/2), where 9 is the 
misorientation angle. This was later verified both in experiments and atomistic simulations [HITT]. In 
general, however, most of the GBs are mixed, containing multiple sets of edge and screw dislocation 
arrays. If the above definition of the geometric coupling factor is generalized to an arbitrary GB, the 
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Figure 2: Schematic of a bicrystal with a planar GB. 


result will be a vector given by 

. PM 

^ V 


(54) 


We will now use this definition to derive an expression for the coupling factor for three special cases: 
(all for a planar GB ^ as shown in Figure [2|) i) Symmetric tilt boundary with finite misorientation. 
Here we will provide an alternate argument to recover the formula obtained earlier by Cahn et al. [3]. 

ii) Twist GB with small misorientation. We show that the coupling factor for such a boundary is zero. 

iii) Mixed GB with small misorientation. 


4.1 Symmetric tilt GB 

Let F be the total deformation gradient of the grains with respect to a fixed reference configuration. 
Gompatibility at the boundary requires {FJ = and |u] = —V^a, where a is an arbitrary vector, 

while rir and Vr, respectively, are the normal vector and the normal velocity of the GB in the reference 
configuration. Without loss of generality, we can assume that F^ = I and = 0. Gonsequently, 
Ur = n and Vr = V m- On the other hand the multiplicative decomposition of F, under the present 
assumption of elastically rigid grains, takes the form F = RF^ [TT], where R is the lattice rotation 
tensor and F^ is the plastic deformation gradient. If we assume the plastic deformation to be isochoric 
(detF'^ = 1), and that F^^ = I, then the above considerations lead to F~ = I+(3®n, where (3 = —Pa 
is a tangential vector (the superscript will be suppressed hereafter). The total Burgers vector B 
of all the GB dislocations cut by a unit vector p lying on the GB plane is given by the Frank-Bilby 
equation [III 

B = {I- R^)p = (/ - FP)p. (55) 

Assuming all the edge dislocations at the GB to glide in a single slip direction, we can write the 
resulting plastic distortion rate as F^{FP)~^ = (^s ® m (cf. Ghapter 106 in [l3]), where C, s, and 
m stand for slip rate, unit slip vector, and unit normal to the slip plane, respectively (s and m are 
mutually perpendicular). With initial values of C and F^ as 0 and I respectively, time integration of 
the evolution equation yields F^ = exp(^(t)s < 8 ) m){FP\t=o) = I + 0 m. If the orientation of grain 

G~ is related to that of G~^ by an anticlockwise rotation of angle 63 about e 3 -axis (see Figure [2]) then 
we can write R = cos ^ 3(61 0 ei + 62 0 62 ) + sin 03 (e 2 0 ei — ei (8 62 ) + 63 (8 63 . Using this in (f55]l i for 
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p = 62 , and recalling that B = |-B|s, we obtain \B\ = 2 sin( 03 / 2 ) and s = — cos( 03 / 2 )ei +sin( 03 / 2 )e 2 
(and hence m = — sin( 03 / 2 )ei — cos( 03 / 2 )e 2 ). Finally, with the help of expressions derived above for 
total and plastic deformation gradients, we can obtain Q = —\B\/{m ■ 62 ) and consequently, 

/3 = 2tan(03/2)e2. (56) 

This expression for the coupling factor was derived earlier by Cahn et al. [3]. For small misorientation 
angle the coupling factor takes a simple form (3 = ^362 (Chapter 14 in [19] 1. 

A general tilt boundary with small misorientation: Restricting ourselves to small misorientation, we 
consider a tilt GB such that R = I + 92 (^ 2 ^) + ^ 3 (^ 3 x), where (ex) represents a skew tensor with 
components given by (ex)jfc = SjikSi (here Sjik is the permutation symbol). In other words, grain G~ 
is obtained by rotating the reference grain anticlockwise about 62 and 63 by small angles 62 and 9 ^, 
respectively. The dislocation density tensor at the GB, dehned as a = (/ — F^){n x) = (/ - R^){nx) 
mm, takes the form 

et = 0261 0 62 + 0361 0 63 . (57) 

This represents two arrays of edge dislocations having line direction along 62 and 63 , and slip direction 
61 , with densities 02 and 03 , respectively. To calculate the geometric coupling factor we exploit linearity 
in extending the above result for a symmetric tilt boundary to the present situation to obtain 

/3 = 0362 - 0263 . (58) 

The coupling factor therefore has contributions from both the arrays of edge dislocation. 

4.2 Twist GB with small misorientation 

We now consider a twist GB with small misorientation such that the grain G~ is rotated by an 
anticlockwise angle 0i, about 61 -axis, with respect to grain G"*". For small angle we can write R = 
J -I- 01 (61 X). As a result 

Q = - 0 l (62 (g) 62 + 63 (g) 63 ), (59) 

which represents two arrays of screw dislocations (both with density 0 i) with line directions parallel 
to 62 and 63 . With this in mind we assume the plastic deformation gradient as = I + ^ + 

S 2 (g m 2 ), where we have considered two mutually-orthogonal slip systems with equal slip magnitude, 
such that |si| = |s 2 | = |mi| = |m 2 | = 1 and si • S 2 = ■ mi = S 2 ■ m 2 = mi • m 2 = 0. On the 

other hand, the total deformation gradient is of the form considered above, i.e. F = 7 -|- /3 (g 61 . With 

the assumption of small deformation and small misorientation, the multiplicative decomposition of the 
deformation gradient becomes an additive decomposition so as to yield the following for the case at 
hand: 

/3 g) 61 = 0i(6ix) + C(si g) mi -b S 2 g> m 2 ). (60) 

Projecting this onto 61 , 62 , and 63 we obtain 

/3 = C(si(mi • 61)-bS2(m2 • 61)), ( 61 ) 

0 = 0163 -b C (si(mi • 62) + S2{m2 • 62)), and ( 62 ) 

0 = -0162 + C (si(mi • 63) -b S2(m2 • 63)), ( 63 ) 
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respectively. After some manipulations, the latter two equations yield si ■ ei = S2 ■ ei = 0, mi = S2, 
and m 2 = si- These results, when substituted into (j61|) . immediately furnish /3 = 0. We have therefore 
shown that the geometric coupling factor for a twist boundary (with small misorientation) is zero, hence 
confirming the qualitative arguments provided in [3l[5l[Tn] . 

4.3 A cubic grain embedded in a large grain 

As an application of the results obtained in the previous two subsections we now consider an example 
where a cubic grain (whose edges are aligned with directions ei, 62 , and 63 ) is embedded inside another 
grain such that the (inhnitesimal) misorientation between them is given by 

R = I + 9i{eix) + 02(e2x) + 6'3(e3x). (64) 

The surface dislocation density tensor for the GB with normal ei can be calculated as 

Q = -9i{e2 ( 8 ) 62 + 63 (g) 63 ) + 0261 (g) 62 + 6*361 (g 63 , (65) 

which is the sum of densities given in (I57|) and (I59p ; the GB is of a mixed type consisting of two mutually 
perpendicular sets of edge dislocations (with densities 62 and 6 ^) and two mutually perpendicular sets of 
screw dislocations (both with densities 61 ). The dislocation content at other boundaries can be obtained 
in a similar manner. In evaluating the geometric coupling factor associated with the boundary with 
normal 61 , we exploit linearity in our arguments (due to small misorientation) to combine the results 
obtained above for tilt and twist boundaries to write 

/ 3 ei = ^362 - 6*263. (66) 

We can similarly calculate the coupling factors for the GBs with normal 62 and 63 as 

= “ 6*3 61 + 6*1 63 , and / 3 e 3 = 6*2 61 - 6*1 62 , (67) 

respectively. It is easily verihable that /3_ej = —etc, where /3_e^ represents the coupling factor 
associated with the face with normal — 63 . 


5 Kinetic relations 


The governing equations for coupled GB dynamics with junctions can be derived starting from in¬ 
equalities (1471) and (1451) by first identifying various dissipative fluxes, and the associated driving forces, 
and then assuming linear kinetics. Towards this end we consider three crystalline arrangements: (i) 
bicrystal-I consisting of two rectangular grains joined at a mixed planar GB (as in Figure [2]) and sub¬ 
jected to shear stress; (ii) bicrystal-II with a spherical grain embedded inside a larger grain (see Figure 


3(a)); and (hi) tricrystal where a grain Gi is embedded inside a large bicrystal made of two rectangular 


grains G 2 and G 3 (see Figure [3(b)[). 
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5.1 Bicrystal-I 

The first bicrystalline arrangement is as shown in Figure[2]such that the traction on the outer boundaries 
perpendicular to ei is Te 2 and the traction on the outer boundaries perpendicular to 62 is rei. We 
also assume that the grains are rigid, free of defects, and contain negligible stored energy. Moreover, 
we neglect all kinds of atomic diffusion. The GB 6 ^ is considered to be of a mixed type, with both tilt 
and twist components, where the misorientation is given by 

© = 7 +6li(eix) + 6l3(e3x), (68) 

for small 61 and ^ 3 ; as discussed previously, 9i and 03 determine the twist and the tilt characteristic, 
respectively, of the GB. Whereas the deformation of grains in the wake of a moving GB, under the 
external loading considered here, is simple shear for a tilt GB, it is more complicated if the GB is of 
mixed type |lf)] . The array of edge dislocations are driven by the Peach-Koehler force to move the GB 
in normal direction while translating the grains parallel to the GB. On the other hand, the simultaneous 
movement of two perpendicular sets of screw dislocation arrays results into a relative rotation of the 
adjacent grains about the GB normal. The GB motion, the relative tangential translation, and the 
grain rotation are in general all coupled to each other. 

The state of stress throughout the bicrystal is taken as cr = r(e 2 <Xiei + ei( 8 )e 2 ); this clearly satisfies 
both the equilibrium equations and the traction boundary conditions. Based on the experimental 
observations in unmn, we assume the tilt angle to remain hxed while allowing the twist angle to 
evolve owing to the relative rotation between the grains. The axis of rotation is taken to coincide with 
63 . Without loss of generality, the grain can be assumed to remain stationary, i.e. v~^ = 0, and 
G~ moving with a velocity 

v~ = 0iei X X + Ce2, (69) 

where G is the translational velocity of G~ in the direction of 62 and x = xiei + 0:262 + 0:363 is the 
position vector. Observing that lunl = 0, the dissipation inequality (ITFl) reduces to 

0i/0 + G/c>O, (70) 

where fg = {tx^ — d'y/dOi) and fc = —r. Considering Onsager’s reciprocity theorem [18], we can 
obtain the following pair of coupled kinetic relations | 2 |: 

01 = Safe + BSdfc and C = 861 + Cf„ (71) 

where <5^ > 0 is the sliding coefficient due to the relative rotational motion (caused by the intrinsic 
screw dislocation glide along =5^), T > 0 is the sliding coefficient for the relative translational motion 
between the grains, and B denotes a coupling between grain rotation and translation. In response to 
grain translation, the edge dislocation array will cause simultaneous GB migration |3|, such that 

V = -C/h. (72) 

where (32 = 03 is the geometric coupling factor as calculated in the previous section. Note that the 
geometric coupling exists only with respect to the translational velocity, as the rotational part in (16911 


15 




Figure 3: (a) Schematic of bicrystal-II. (b) Schematic of tricrystal. 


amount to pure sliding. According to (j72p it is the direction of C which decides whether the GB will 
move upwards or downwards. 

To summarize, we have a coupled system of equations, given by m and (172]) . which should be 
solved to evaluate the position of the grains and the GB, as well as the misorientation, at any given 
time instance during the dynamical process. It should be noted that, in a more complicated situation 
when the driving forces are functions of xi and X 2 , an initially planar GB will not necessarily remain 
planar (cf. |10] 1 and GB diffusion will be required to prevent void-formation/interpenetration at the 
GB. 


5.2 Bicrystal-II 

As a second example, we consider a bicrystal with a spherical grain G~ (of radius R) embedded within 
a much larger grain G~^, as shown in Figure 3(a), with misorientation between the grains given by 

© = J + 0i(eix) + 02(e2x) + e^{e^x), (73) 


i.e. grain G~ has been obtained by rotating it from G"*" by small angles 62 , and 63 about ei, 62 , and 
63 , respectively, where the orthonormal basis vectors { 61 , 62 , 63 } form a coordinate frame with origin 
at the center of the sphere. The GB ^ hence is of a mixed type. Without loss of generality we let 
the outer grain G^ to remain fixed and allow the inner grain G~ to rotate (without translating). We 
assume that the external stress is absent, the grains are rigid and free of defects, the free energy of the 
grains is vanishing, and the volumetric diffusion is absent. Let us consider an orthonormal spherical 
basis {eR,e^,e^} with origin at the center of the embedded grain (^ is the polar angle and (p is the 
azimuthal angle). The GB normal n points into G'*', hence n = en. The angular velocity of G~ (axial 
vector of 0 ) is given by 

m = 01 61 + 02 £2 + ^3 63 - (74) 

Neglecting rigid body translation, the velocity of the inner grain can be written as v~ = w x x, where 
X is the position vector. Recalling that v~^ = 0 we obtain 


bnl = Vn = 0 and Vt = P(^»+ -v ) = -v , 


(75) 
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where Vt represents the relative tangential velocity of G~ with respect to G~^ at the GB. 

The spherical GB rotates and shrinks without any shape change. This does not require any shape 
accommodation mechanism such as GB diffusion. We will therefore consider h = 0. Assuming isotropic 
GB energy and defining 6 = 9iei + 02^2 + 6 * 3^3 (the axial vector of © — /) the dissipation inequality, 
given by (H7D . reduces to 

Vfn + w ■ > 0 , where (76) 

fn = IK and ft = -de'j. (77) 

Invoking the Onsager reciprocity theorem m, we can postulate the following linear kinetic relations 
based on dZSD 12]: 

V = Mfn + M^-ftS^^d (78) 

w = pV + S~ft, (79) 

where Ad > 0 is the GB mobility, f3 is the coupling factor between rotational speed and normal GB 
velocity, and S is the symmetric positive semi-definite sliding coefficient. Both the viscous effect and 
the twist characteristic of the GB are expected to contribute to the net sliding. To understand the 
physical meaning of the coefficients /3 and S, we take a cross-product of ()79p with x to obtain 

Vt = f3V + Sft, where (80) 

(3 = xxP, Sft = xx Sft, and ft = {l/R)eR x ff (81) 

Equation (| 8 ip i relates (3 to the the geometric coupling factor f3 introduced in the previous section, while 
m 2 relates S to the sliding coefficient S; the relation dHIjls is motivated from the 2D counterpart of 
the present discussion [2]. Next, we specialize these kinetic equations for a spherical GB, under various 
additional assumptions, and present analytical solutions wherever possible. When reduced to a 2D 
setting, the derived relations will be identical to those obtained for a circular GB in [S] and [2]. 

GB migration: When both geometric coupling and GB sliding are negligible, kinetic equations (I78p and 
()79p simplify to the well-known equations for curvature driven GB migration: V = A4jk and w = 0. 
Using V = R and k = —1/R (for a spherical GB) in the governing equations we obtain the following 
solutions: 

R{t) = Rq — 2M.'yt and 6{t) = 6q, (82) 

where Rq is the initial GB radius and Oq is the initial misorientation. Both Ad and 7 have been treated 
as constants. 

Coupled motion without GB sliding: At temperatures far below the melting point, GB viscous sliding 
becomes negligible and geometric coupling plays the dominant role in grain rotation (3|. The kinetic 
relations (fTsp and ([7^ can then be assumed to take the form 

V = Mjk + Mfi ■ ft w = PV, (83) 

respectively. For a spherical GB (with V = R and (3 = {& — (ISiP i and a cross product of (fMP o 

with X = Rsr furnishes 

R 

wxeR = -—9xeR, (84) 
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which on time integration gives 


R{t)e{t) = RqOo 


(85) 


for all ^ and cj). Note that we have used the same geometric coupling factor associated with a planar GB 
for analysing the kinetics of a curved GB; this assumption is motivated from the 2D atomistic studies 
where the geometric coupling factor was observed to be same for planar and curved GBs (cf. [ 22 ] and 
the references therein). Before integrating (j83|) i we consider the following isotropic energy proposed 
by Read (see Section 12.7 in [H]): 

7 = 7 o| 0 |(Ac-ln| 0 |), (86) 

where 70 is a constant depending on the material properties, the Burgers vector, and the spacing be¬ 
tween the dislocations; Ac is a constant which depends on the energy of atomic misfit at the dislocation 
core. Differentiating (I 86 p with respect to 6 we obtain the torque 

=-007 =- 7 o|^(Ac - 1 - ln|0|). (87) 

Substituting this into (l83]) i . and restricting it to a spherical GB (for which V = R, k, = —1/R, and 
(3 = —6/R), we obtain R = —Al7o|0|/7?. Gombining it with ([851) and subsequently integrating yields 

R{t) = (Rq - 3MjoRo\6o\t)^/^ and e{t) = Ro{RI - 3M-foRo\Oo\t)~^/^Go, ( 88 ) 

where A4 has been treated as a constant. 

Coupled motion without geometric coupling: At temperatures close to the melting point, viscous sliding 
at the GB dominates over the geometric coupling to govern grain rotation |3]. The kinetic relations 
(f78]l and ([7^ in such case reduce down to 

V = Mjk, and G = Sf^, (89) 

respectively. Solving these equations analytically for R and G is challenging; we will attempt to derive 
a useful implicit relation between them. If we assume S and S to be of the form S = SP and S = SI, 
respectively, then (f 8 T]l 9 requires S = S/R^. Gonsequently, for a spherical GB, equations (l89l) can be 
manipulated to obtain 

0 ^ 5 G{Ac-l-ln\G\) 

R MR |6»|2(yl^-ln| 6 »|) ’ ^ ’ 

where we have also used (| 86 l) and (|87l) . This is a non-linear equation in G. Taking a dot product with G 
on both sides of the equation, and then integrating the result, we obtain the following implicit relation 
between R and G: 

R{G) = Roexp {|0|2 - \Go\^ + {I, u^) - R(l,u)]}^ , (91) 

where u = 2{Ac — 1 — ln|0|), uq = 2(Ac — 1 — ln|0o|)) and E{n,y) = yi^-n^n du (the exponential 
integral); M and S have been considered to be constants. 
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Fully coupled motion: We finally consider the situation when both geometrical coupling and GB sliding 
will contribute comparably to the evolution of grain rotation. Assuming S to be invertible, and 
eliminating between (f78ll and d79]l . we derive for a spherical GB (/3 = —6/R, V = R, and k = —1/R) 


R = - ^(Mj + MS ^e -w 

R^ + Me-s 

An expression for w can be obtained by substituting (I78]l in ([7^ : 


(92) 


M .^6 


~ M 


0 ft 


(93) 


5.3 Tricrystal 

We consider a tricrystal, as shown in Figure [3(b)[ where grain Gi is embedded inside a larger bicrys¬ 
tal made of two rectangular grains G 2 and G 3 . The tricrystal is subjected to external stress. The 
configuration has three GBs with unit normals denoted by rij {i = 1,2,3), and a closed junction 
curve J. The normals are chosen such that both rii and n2 point into Gi whereas 723 points into G2. 
In deriving the kinetic laws we assume the following: (i) the grains are rigid, defect free, and have a 
vanishing stored energy; (ii) volumetric diffusion is ignored; (iii) the shape accommodation required for 
preventing void-formation/interpenetration at various GBs is accomplished by allowing for diffusion 
along the GBs; (iv) diffusion along the junction is negligible; and (v) the magnitude of applied stresses 
are small enough so that elastic and plastic deformation of the grains can be neglected. Under the 
combined effects of GB capillary force and the applied stress field the GBs will migrate, grain Gi will 
rotate and translate (as a rigid body), grains G 2 and G 3 will translate rigidly relative to each other, 
and the junction J will move in space, all coupled to each other. 

In the absence of intra-granular defects the orientation field associated with grain Gi, denoted by 
Ri, will be homogeneous throughout the grain. We define the misorientation tensor for the respective 
GBs as ©1 = R[R2, ©2 = RJR3, and ©3 = R2R3. We assume grains G2 and G3 to be non-rotating, 
thereby fixing their orientations once for all. As a result 

0i©f = ©2©!^ = -RiRj. (94) 

On the other hand, the velocities of rigidly deforming grains can be written as 

vi = w X X + Cl, V 2 = C 2 , V 3 = 0, (95) 

where w is the angular velocity of Gi (the axial vector of ), x is the position vector, and Ci is the 

rigid translation of Gi (grain G 3 has been assumed to remain fixed). We define the relative translation 
velocity of adjacent grains at the respective GBs as Ci = Ci — C 2 , C 2 = Ci, and C 3 = C 2 . Before 
we substitute these in the mass balance relations, we would additionally assume that the translational 
velocity of the embedded grain to be much smaller than its rotational velocity. We can justify this 
on the basis of the atomistic simulation results which do not show any significant translation in the 
absence of external stress |22] . The small amplitude of external stress considered here would therefore 
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cause only small relative translation. Keeping this in mind and using (I95p in (I24p . along with negligible 
diffusional fluxes in the grain, we derive 

div^ ha = —px X Ua ■ w for a = 1, 2, and 

div'^ /13 = —pCs ■ ns. (96) 

These relations can be integrated and then combined with (j25p (where j = 0, hj = 0, and p, v are 
non-singular at the junction) to obtain 

ha = AaW and h^ = A 3 C 3 , (97) 


where Aj is a second order tensor which depends on the geometry of =5^, (see [T] for a similar calculation 
for a 2D tricrystal). We also relate these fluxes to the chemical potential by assuming the Tick’s law 
for superficial diffusion [T3] 

hi = —DiV^p for z = 1, 2,3, (98) 

where Di is the (symmetric and tangential) diffusivity tensor along .5^. 

For the present case, the dissipation inequality in the bulk (j46l) is trivially satisfied. The dissipation 
inequalities at the GBs, given by (j47p . are however non-trivial and will be used in the following to 
derive the kinetic equations. Substituting (IMp . ([95]), (IWP . and ([MP into (H7p . and neglecting the terms 
of the order of |CaP and |i?aCa|, we can reduce the inequality to the form (no summation for repeated 
index i) 

Vifi + w ■ Gi + Ci ■ Hi > 0 on y’i, (99) 

where 

Hi = arii, Ga = XX {aUa -h ppria) -h A^Da^AaW + 2Ta, and ^3 = 0. (100) 

In the above relations, is the Moore-Penrose pseudoinverse of Di which satisfies D^^Di = 

DiD-^ = Pirn, and Tq is the axial vector of skew part of {d&^'ya)®a ■ Considering the On- 
sager’s reciprocity theorem [18], we use (IMP to postulate linear kinetic relations associated with S^a 
(a = 1, 2). In writing them we assume the relative translational velocities Ca to be decoupled from GB 
migration and the rotation rate of grain Gi; a theory without this assumption can easily be constructed 
along similar lines. The kinetic relations are taken as 


Va=M[^^fa + Mi^^ -Ga, ( 101 ) 

W = M^^^fa + M^3^Ga, ( 102 ) 

Ca = CaHa, (103) 

where Ca are various kinetic coefficients, to be discussed next. Assuming 

to be non-vanishing we eliminate fa from ( 11021 ) using (IIUII) to rewrite w as 


W = ^aVa + SaGa, (104) 

where M.a = \ Pa = /Ma, and Sa = — A4aPa ® Pa- Th® coefficients Ada, Pa-! and Sa 

have the same physical interpretation as described in Section 15.21 Substituting (|10ip , (|103p , and (I104p 
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back into the dissipation inequality (|M1) we derive the following restrictions: Ma > 0; Sa and Ca are 
symmetric positive semi-definite. In terms of these new kinetic coefficients (jlOip takes the form 

I4 = >fa(/a + ^a'C^a). (105) 

Furthermore, when sliding is active and Sa is invertible, we can eliminate Qa from (11051) with the help 
of (llOdp to obtain 

K = - (fa + Pa • . (106) 

l+MaPa-Sa ^ 

Using (llOOp and (jl05p in (I104p . and rearranging the resulting expression, we can derive 

(^I - ZaA^D^^Aa^ w = MafaPa + ^a{Xa X Crria 2Ta + P^Xa X n^), (107) 

where Za = Sa + MaPa ® Pa- For an invertible Sa we multiply both sides of (I107p by Z^ ^ and 
integrate the result over and 5^2 for a = 1 and a = 2, respectively. We add the two expressions to 
obtain the following expression for w. 

'W = ('^ [ i^a^ - AaDa^Aa)da\ '^f (^Mafa^a^Pa + XaXO-na + 2 Ta'jda, (108) 

where we have used fy PPXa x Ua da = p x x Vpdv = 0 (the hrst equality follows from the 

divergence theorem, whereas the second equality holds due to vanishing diffusional flux in the grain). 
To summarize the results obtained so far, we have the governing equations for the motion of and 
^2 in (llOOp . and the governing equation for rotation of Gi in (llOSp . 

The kinetic relations for ^3 can be derived similarly. In doing so, however, we allow for translational 
velocity to couple with the normal motion. Starting with (|99p . and assuming linear kinetics, we obtain 

C 3 = Ps^s + -^ 3 ^ 3 , (109) 

V3 = M3(f3+p3-n3), ( 110 ) 

where Ms > 0 , P3, and £.3 (positive semi-dehnite) are the mobility, geometric coupling factor, and 
sliding coefficient, respectively, for ^ 3 . Replacing V 3 from (IllOl) in (11091) the expression for translational 
velocity can be rewritten as 

C3 = MsfsPs + ZsHs, (111) 

where Z3 = C3 + M3P3®P3. Since the translational velocities are homogeneous, and the outer grains 
have been assumed to be much larger than the embedded grain, we integrate (lllip over ^3 to write 

C3= [ {M3hP3 + Z3'H3)da. (112) 

On the other hand, the governing equation for average translation velocity of the embedded grain can 
be obtained by first integrating (11031) for a = 1 and 2, respectively, and then adding them to obtain 

= »rea(.4uJ1,) + g CM.dPj , (113) 
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where C 3 is given by (I112p . In (jllOp we have the governing equation for the normal motion of = 5^3 and 
in (lll.ip for the translation of the embedded grain. 

Finally, we derive the kinetic relations which govern junction dynamics. For negligible diffusion 
along the junction curve, the dissipation inequality ([5^ simplihes to > 0. Assuming linear 

kinetics we postulate that 

Qp = M-jTj, (114) 

where Ad .7 is the positive semi-definite junction mobility tensor. An analogous treatment in a 2D 
setting can be seen in [ 8 ] (see also m)- The junction force iFj given by (f53P is a function of the 
unknown local orientations of the adjacent GBs which, for a non-splitting junction, can be calculated 
using the compatibility conditions V) = Qp • (see m for a detailed calculation in 2D). 

To conclude, the complete set of kinetic equations governing the coupled GB motion in the tricrys¬ 
talline arrangement includes ()106p for the motion of and ^5^2; dHOP for the motion of = 5 ^ 3 , (llOSp 
for the rotation of the embedded grain Gi (the outer grains are non-rotating), (I113p and ()112p for the 
translation of grains Gi and G 2 , respectively (whereas grain G 3 is stationary), and (|114p (in association 
with the compatibility condition) for the motion of the junction curve. 

6 Concluding remarks 

We have presented a thermodynamically consistent 3D study of coupled GB motion in the presence 
of junctions, hitherto restricted to 2D crystalline materials. Towards this end we introduced a novel 
continuum mechanics based theory of irreversible dynamics of incoherent interfaces with junctions, 
which allows for diffusion in the bulk, on the interface, and along the junction curve. The various 
local dissipation inequalities derived therein were used to motivate kinetic relations for the coupled 
GB motion in two bicrystals and one tricrystal. These relations were solved analytically whenever 
it was possible to do so, but were otherwise left in a form amenable to numerical computations. In 
any case, the results clearly demonstrated the effect of coupling on the grain dynamics. Gonsider for 
instance the shrinking of an isolated grain, embedded within a larger grain, under the action of capillary. 
Without coupling, the embedded grain can disappear only by shrinking to a vanishing size. However 
with coupling, the grain can disappear by aligning its orientation with the outer grain even before it 
has shrunk significantly [2]. The proposed kinetic relations also emphasize the coupling of junction 
dynamics with both grain and GB motion. Depending on the junction mobility, grain dynamics can 
experience a substantial drag compared to the case with no junctions [T]. 

The present work can form a basis for research in several future directions. The theory of incoherent 
interfaces, which includes junctions and diffusion, is in fact applicable to a more general situation where 
the grains are allowed to deform plastically. The resulting framework would be useful for phenomena 
which involves coupling of plastically deforming bulk with moving incoherent interfaces and junctions. 
Secondly, the kinetic relations derived for the tricrystal can be used to study the coupled motion in 
polycrystalline materials containing large number of grains and junctions. Of course, this will demand 
signihcant computational effort and hence efficient numerical algorithms. A related direction of work 
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would be to develop numerical techniques (such as level set methods) for solving equations of coupled 
motion of an embedded grain with anisotropic constitutive properties. Thirdly, the computation of 
the vectorial coupling factor, which has been restricted here to small angle GBs, should be extended 
to large angle GBs. Finally, 3D atomistic simulations will be required to clarify the nature of various 
kinetic coefficients (including the coupling factor), in particular regarding their dependence on three 
misorientation angles and two orientation angles of the GB. 
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